Astron. Nachr. / AN , No. , 1 -HO / DOI please set DOI! 



Towards an Accurate Model for the Antennae Galaxies 

SJ. Karl^ % T. Naab\ P.H. Johansson\ Ch. Theis^, and CM. Boily^ 

^ Universitats-Sternwarte Munchen, Scheinerstr. 1, D-81679 Miinchen, Germany 

^ Institut fiir Astronomie der Universitat Wien, Turkenschanzstr. 17, A- 1 180 Vienna, Austria 

^ Observatoire astronomique, Universite de Strasbourg and CNRS UMR 7550, 11 rue de I'Universite, F-67000 Stras- 
bourg, France 



00 

o 
o 

<N 
CO 

a^ 

<N 

Oh- 

o 



> 

O 
(N 
O 

in 

On 
O 
00 
O 



The dates of receipt and acceptance should be inserted later 

Key words galaxies: evolution - galaxies: interaction - galaxies: individual (NGC 4038, NGC 4039) - 
methods: numerical 

In the framework of hierarchical structure formation ellipticals can form from merging of smaller disk galaxies. The 
nearby interacting 'Antennae' galaxy pair (NGC 4038/39) is one of the best-studied local systems of merging spirals, thus 
presenting us with an ideal laboratory for the study of galaxy evolution models. The Antennae are believed to be in a state 
prior to their final encounter with rapid subsequent merging, which puts them in the first position of the Toomre (1977) 
merger sequence. Here we present first numerical high-resolution, self-consistent, smoothed particle hydrodynamics (SPH) 
simulations of the Antennae system, including star formation and stellar feedback, and compare our results to VLA HI 
observations by Hibbard et al. (2001). We are able to obtain a close, but not yet perfect match to the observed morphology 
and kinematics of the system. 
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1 Introduction 

The 'Antennae' galaxies (NGC 4038/39, Arp 244, VV245) 
are a well-known, archetypal example of a spiral-spiral 
galaxy merger (see Fig.[T]). They are nicknamed after their 
spectacular appearance with a prominent pair of tidal tails, 
which was already noted in early optical images (e.g. Dun- 
can 1923, Schweizer 1978). At present there is a huge 
amount of data collected for the Antennae both from 
ground- and space-based telescopes, covering a wide range 
of wavelength regimes. Recent examples include: optical 
HST WFPC data (Whitmore et al. 1999), near-IR WIRC 
(Brandl et al. 2005) and mid-IR IRAC (Wang et al. 2004) 
imaging, and Chandra ACIS-S X-ray observations (Baldi et 
al. 2006). The elongated tails are a clear sign of tidal in- 
teraction of nearly equal-mass galaxies (Toomre & Toomre 
1972). They form kinematically by gravitational tides ex- 
erted during the interaction process and their morphology 
and velocity fields give strong hints on the encounter geom- 
etry and history. This makes the Antennae galaxies a key 
system for investigating galaxy interactions and their asso- 
ciated physical phenomena. Starting from the first numer- 
ical model by Toomre & Toomre (1972), where they used 
the restricted N-body method to model the evolution of the 
disks, there have been many attempts to obtain numerical 
'look-alikes' for the Antennae galaxies. Barnes (1988) was 
the first to use self-consistent N-body models with multiple 
components consisting of a bulge, a disk, and dark halo with 
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a mass ratio of 1:3:16 and a total mass of 2.75 • 10^^ M©. 
Dubinski, Mihos & Hernquist (1996) used the extent of the 
tidal features in the Antennae to probe the amount of dark 
matter in these galaxies. Mihos, Bothun & Richstone (1993) 
were the first to include gas dynamics and star formation 
into their dynamical model of the Antennae galaxies. For 
the latest review of the Antennae models we would also like 
to refer to the results presented by Hibbard (2003). How- 
ever, most of the Antennae models are still based on the 
orbital parameters given in Toomre & Toomre (1972). Here 
we want to extend the probed parameter space in order to 
find orbital parameters which allow for a closer fit to the 
observed morphological and kinematical data of the An- 
tennae. We use medium-resolution VLA HI mappings of 
the Antennae (-- 20", = 5.21km s"^) by Hibbard et 
al. (2001) for comparison with our models. In this case it is 
advantageous to use the cold atomic gas as a sensitive tracer 
of the large-scale dynamics of the system as it is unlikely 
to be disturbed by star formation. Recently, there has been 
a debate about the exact distance to the Antennae, ranging 
from a modest 13.3 ±1.0 Mpc (Saviane et al. 2008), based 
on photometry of the tip of the red giant branch, to 22 ± 3 
Mpc (Schweizer et al. 2008), based on observations of a su- 
pernova of type la. Sometimes even higher and lower values 
have been quoted (e.g. Zezas & Fabbiano 2002, Rubin, Ford 
& D'Odorico 1970). A widely-used, intermediate distance 
to the Antennae is I) = 19.2 Mpc, which is derived from 
the systemic recession velocity relative to the Local Group 
assuming a Hubble constant of Hq = 75 km Mpc~^ 
(Whitmore et al. 1999). In this paper the distance will be 
adopted as part of the model matching process. 
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Fig. 1 NGC 4038 (south) and NGC 4039 (north). Distri- 
bution of HI and optical light. HI data are shown in blue, 
together with a combined B + V + R optical image in white 
and green. Taken from Hibbard et al. (2001). North is point- 
ing upwards. 
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Fig. 2 'Best match' to the Antennae at t = 0.60 Gyr in 
the plane-of the sky. Disk and bulge particles are shown in 
green, gas in blue, and stars, which have formed during the 
simulation, in red. Note that for the old stellar component 
(green) we plot only every fifth particle. 

2 Numerical setup 

For our high-resolution model of the Antennae galaxies we 
set up equilibrium galaxy models following Springel, di 
Matteo & Hernquist (2005). Each model consists of a NFW 
halo (Navarro, Frenk & White 1997) which is then con- 
verted to a Hernquist (1990) profile dark matter halo. Em- 
bedded in the dark halo is a stellar Hernquist (1990) bulge 
component, and an additional exponential stellar and 
gaseous disk. The ratio of luminous to dark matter, or, the 
baryon mass fraction, is fb = Mbary/^DM = 1/4, where 
NGC 4038 is modeled as an 'Sc-type galaxy' with a bulge- 
to-disk ratio of B/D = 0.2 and NGC 4039 as an 'Sb-type 
galaxy' with a B/D = 0.4. Each galaxy has a total mass of 
Mtot = 2.77- lO^^M0.Asa starting point for our models 
we adopted the baryon mass fraction and the total masses 
from Barnes (1988). The initial gas fraction in the disks was 
chosen to be 20%, which amounts to M^^^^ 9.2 • 10^ 
and Mg^g^^ 7.9 • 10^ M©, respectively. Further param- 
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Table 1 Galaxy model parameters 



Modeled Property 


NGC 4038 


NGC 4039 


Disk scale length 


3.31 kpc 


3.29 kpc 


Disk scale height 


0.66 kpc 


0.66 kpc 


Bulge scale length 


0.66 kpc 


0.66 kpc 


Maximum rotational velocity 


169 km/s 


170 km/s 



eters for the two model galaxies are summarized in Tab. 
[U All simulations are run on the local ALTIX SGI super- 
computer using the fully parallel smoothed particle hydro- 
dynamics (SPH) (see e.g. Monaghan 1992) code Gadget2 
(Springel 2005). The simulations include a prescription of 
radiative cooling for primordial hydrogen and helium (Katz, 
Weinberg & Hernquist 1996). We also model star formation 
and supernova feedback, following the sub-grid multiphase 
prescription as described by Springel & Hernquist (2003), 
but exclude supernovae-driven galactic winds. Initially the 
galaxies move on elliptical Keplerian orbits (e ^ 0.8) with 
a pericentric separation rp = 20 kpc and an initial sepa- 
ration of one virial radius rinit = 106 kpc. We ran a set 
of low-resolution simulations, going through iterative cy- 
cles, while varying the scaling parameters, the orientation 
of the disks and the viewing angle, until we obtained a good 
match with the observational data. However, this resulted in 
a rather high adopted distance to the Antennae of D = 32.6 
Mpc if we choose to apply no spatial scaling. This merger 
geometry was then re- simulated at high-resolution includ- 
ing star formation and radiative cooling. The total number 
of particles was A^tot = 1,600,000, of which NGC 4038 
and NGC4039 contributed 105,000 and 180,000 bulge par- 
ticles, 420,000 and 360,000 disk particles, and 105,000 and 
90,000 gas particles, respectively. Furthermore, each galaxy 
consisted of A^dm = 170,000 dark matter halo particles. 
The particle numbers were chosen in this particular way 
for dynamical reasons, i.e. to ensure that all particles in the 
baryonic component have exactly the same masses in order 
to minimize two-body relaxation effects. The gravitational 
softening parameters for stellar and gaseous particles were 
set to e = 0.020 kpc and for the dark matter halo particles to 
e = 0.083 kpc, and the system was evolved for a total time 
of - 2 Gyr. 

3 Results 

Fig.[3]shows the evolution of the Antennae merger morphol- 
ogy as a function of time. In our simulation, we set the ori- 
gin of the time such that the galaxies pass pericenter at a 
time t = Gyr (upper panel in Fig. O, and have their fi- 
nal encounter at t ~ 0.82 Gyr. Star formation is mostly 
confined to the central regions of the simulated galaxies. 
At the time we stopped our calculations, t lA Gyr, the 
galaxies still show elongated tidal tails and clear signs of 
a tidally disturbed morphology at their centers. They have 
not yet transformed into an elliptical-like merger remnant 
(see, e.g. Naab & Burkert 2003, Naab, Jesseit & Burkert 
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Fig. 4 Position- velocity diagram of the Antennae model at time of 'best fit' : X-Z (plane of the sky, left), X-Vy (middle), 
and Vy-Z (right). Only gas particles for both simulated galaxies NGC 4038 (red) and NGC 4039 (blue) are shown. HI data 
(Hibbard et al. 2001) are overlaid in yellow. 




Fig. 3 Dynamical evolution in the NGC 4038/39 simula- 
tions. We show a sequence of snapshots at three equally- 
spaced times between pericenter (t = 0.0 Gyr) and best 
match (t = 0.82 Gyr). Colors are coded the same way as in 
Fig. El 



2006, Johansson, Naab & Burkert 2008). We obtain our 
'best fit' at a time shortly before the second encounter, that 
is t ~ 600 Myr after pericenter and t ~ 220 Myr before 
the final merger (see Figs. [2] and lower panel in Fig. O. In 
the position- velocity diagram (Fig. |4]) we see three projec- 
tions in the a) X-Z (plane-of-the-sky), b) X-Vy and c) Vy-Z 
planes. Only gas particles are shown, where blue and red 
particles indicate particles from NGC 4039 and 4038, re- 
spectively. Comparing the model with observations (over- 
laid in yellow), we obtain a good fit with the observed mor- 
phology of the system (left panel in Fig. |4]). However, the 
agreement with the observed velocities of the northern tail 
could be improved. In this region the data show overall 
smaller values (right panel in Fig|4]). This could be achieved 
by slight changes of the orientation of NGC 4039 with re- 
spect to the orbital plane and/or by choosing a slightly dif- 
ferent viewing angle. Also the spatial and kinematical data 
at the bending end of the southern tail still need adapta- 
tion (middle/right panel in Fig. |4]), which may be achieved 
by adopting a flat distribution for the initial gas disks. Our 
model is similar to the model by Hibbard (2003), but in ad- 
dition we included gas dynamics, star formation, and feed- 
back in our simulations. 

In contrast to observations, we do not see enhanced star 
formation in the 'overlap' region at the time of 'best fit' as 
is observed in the Antennae (see, e.g. Wang et al. 2004). 
This may be due to the fact that a significant amount of 
gas has already been consumed by the starburst during the 
first encounter. This will be subject to further investigations. 
Fig. [5] shows the global star formation rate (SFR) of the 
simulated Antennae system. The star formation history in 
the two galaxies is dominated by two major starbursts re- 
sulting from gas funneled to the galactic centers after the 
first (t = 0.0 Gyr, black dot) and second passage (t ^ 
0.82 Gyr, red dot). At the time of best fit (t ^ 0.60 Gyr, 
green dot) we have a SFR of ~ 2.3M0yr~^ which lies 
below the range of observed values. Zhang, Fall & Whit- 
more (2001) report a SFR of 4 < M^/Mq yr"^ < 21 de- 
rived from the total Hq, flux within the disks of the Anten- 
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Fig. 5 Time evolution of the total SFR in the NGC 
4038/4939 model. With filled dots we are indicating the first 
pericenter (t = 0.00 Gyr, black), best match (t = 0.60 Gyr, 
green) , and final merging = 0.82 Gyr, red). 

nae either directly (lower value) or with a correction for ex- 
tinction (upper value). Stanford et al. (1990) give a SFR of 
= 8.8 Mq yr~^ for the total system. Our result is qual- 
itatively consistent with the low values of the total SFR of 
5.4 yr~^ reported in the Antennae simulations by Mi- 
hos, Bothun & Richstone (1993). 

4 Summary & Outlook 

In this work we have presented first steps towards a good nu- 
merical model for the Antennae galaxies. The simulations 
include star formation and stellar feedback using sub-grid 
physics as proposed by Springel & Hernquist (2003). We 
are able to obtain a close match to the observed morphol- 
ogy of the system, while a perfect fit to the kinematical HI 
data is still to be achieved. Our goal is to investigate in de- 
tail the history and distribution of newly formed stars in the 
Antennae galaxies, in particular in the 'overlap' region of 
the two progenitor disks. Close comparison to observations 
then enables us to adjust the star formation prescriptions 
used in numerical simulations and figure out the important 
physical properties controlling star formation in interacting 
galaxies. In this spirit, it is worth examining alternative for- 
mulations for the star formation laws in numerical simula- 
tions, see e.g. Barnes' (2005) shock-induced star formation 
model for the 'Mice' galaxies. There is a number of possi- 
ble processes relevant for the enhanced star formation in the 
Antennae merger, specifically in the 'overlap' region. One 
idea is that collisions of giant molecular clouds (GMCs) 
during the merger process could trigger starbursts (Noguchi 
1991) in interacting galaxies. As the disk filling factor of 
GMCs is rather low (/ < 0.01), Jog & Solomon (1992) 
proposed that a hot ionized, high-pressure gas medium, re- 
sulting from fast collisions between HI clouds, could cause 
a radiative shock compression in the outer layers of pre- 
existing GMCs, thereby inducing bursts of star formation. 
Elmegreen & Efremov (1997) advocate high-pressure re- 
gions, which originate in the large-scale shocks of interact- 
ing galaxies, for being the driver of bound star cluster for- 
mation, while, lately, Renaud et al. (2008) raised the idea 
that compressive tidal modes in galaxy mergers could play 



an important role for the formation of globular clusters. It 
would also be interesting to study whether one can form 
candidates for tidal dwarf galaxies in the tails of the merger 
as is observed in the southern tidal arm of the Antennae 
(Schweizer 1978, Mirabel, Dottori & Lutz 1992, see also 
Wetzstein, Naab & Burkert 2007). 
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